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ABSTRACT 

Supernova (SN) explosions enrich the intra-cluster medium (ICM) both by creating and dispersing 
metals. We introduce a method to measure the number of SNe and relative contribution of Type la 
supernovae (SNe la) and core-collapse supernovae (SNe cc) by directly fitting X-ray spectral observa- 
tions. The method has been implemented as an XSPEC model called snapec. snapec utilizes a single 
temperature thermal plasma code (apec) to model the spectral emission based on metal abundances 
calculated using the latest SN yields from SN la and SN cc explosion models. This approach provides 
a self-consistent single set of uncertainties on the total number of SN explosions and relative fraction 
of SN types in the ICM over the cluster lifetime by directly allowing these parameters to be deter- 
mined by SN yields provided by simulations. We apply our approach to the XMM-Newton European 
Photon Imaging Camera (EPIC), Reflection Grating Spectrometer ( RGS), and 200 ks simulated 
Astro-H observations of a cooling flow cluster, A3112. We find that various sets of SN yields present 
in the literature produce an acceptable fit to the EPIC and RGS spectra of A3112. We infer that 
30.3%±5.4% to 37.1%±7.1% of the total SN explosions are SNe la, and the total number of SN explo- 
sions required to create the observed metals is in the range of (1.06 ± 0.34) x 10 9 to (1.28 ± 0.43) x 10 9 , 
from snapec fits to RGS spectra. These values may be compared to the enrichment expected based 
on well-established empirically-measured SN rates per star formed. The proportions of SNe la and 
SNe cc inferred to have enriched the ICM in the inner 52 kpc of A3112 is consistent with these specific 
rates, if one applies a correction for the metals locked up in stars. At the same time, the inferred level 
of SN enrichment corresponds to a star-to-gas mass ratio that is several times greater than the 10% 
estimated globally for clusters in the A3112 mass range. 

Subject headings: galaxies: clusters: intracluster medium, galaxies: individual (A3112) nucleosyn- 
thesis, abundances supernovae: general X-rays: galaxies: clusters 



1. INTRODUCTION 
The use of Type la supernovae (SNe la) as standardiz- 



able candles in the discovery of dark energy ( Riess et al 



1998 Perlmutter et al. 1999) installed the nature of SIN la 



progenitors and the physics of the explosion as one of the 
prime problems in astrophysics. The evolution of the SN 
la rate is one o f the promising methods for unveiling SN 



la progenitor s ( Ruiz-Lapuente & Canal 1998 Yungelson 



& Livio|2000|). Measuring the evolution of the supernova 



(SJN) rate is particularly important at higher redshifts 



et al. 


2002 


Gal 


axy c 



Barbary et aL|2012| ) 



usters represent the largest scales of orga- 
nized matter in the universe, making them unique lab- 
oratories for chemical enrichment of the universe from 
all possible sources including SN la. The intra-cluster 
medium (ICM) emits X-rays due to the highly ionized 
gas which h as been heated by in fall from the intergalac- 
tic medium ( Gunn fc Gott|[l972 ). The metals produced, 



during all stages of cluster formation and evolution, b 
stars and gal axies enter the ICM via SN e xplosions ([De 



Young| [1978D and strong galactic winds (Mathews 
Baker||197ip . They remain in the ICM due to cluster's 
deep potential well. Since X-ray spectroscopy yields ac- 
curate measurements of metal abundances in the ICM, 
the large reservoir of metals in clusters of galaxies pro- 
vides a unique way to probe SN rates and therefore SN 



la progenitor models on a universal scale. 

The launch of ASCA triggered studies of the rela- 
tive contributions of SNe la and core-collapse supernovae 
(SNe cc) to ICM metal enrichment. These studies used 
specific elemental abundance ratios (e.g., nickel to iron) 
measured in clusters to distinguish between SN la models 
flBaumgartner et al.|[2005l |Dupke k, White||2000l |Dupke| 
fc Arn aud 2001} |Mushotzky fc Loewenstei n 1997). With 
the launch of observatories with better spectral and spa- 
tial resolution, such as XMM-Newton, Chandra, and 
Suzaku, radial profiles of specific elemental abundances 
were used to determine the relative contribution of SNe 
la to the enrichment of cluste r s as a function of position 
e ^Matsushita et al.| ( |2007[ ); |Sato et al.| ( |2007| ); |Million 



et al. 



( 201ip ). Higli^solut ion XMM-Newton Reflection 



Grating Spectrometer ( RGS) measurements of specific 
elemental abundance ratios within the ICM have also 
been used to determine the integral yield of all the SN 
la and SN cc explosions and to distinguish between SN 
models with different levels of pre-enrichment of the pro- 
genitors and with different initial-mass functions (IMF) 
jWgrner et al]|2006b[ |de Plaa et aT1|2007| [Grange et al. 



2011) 



In this paper we introduce a new approach to deter- 
mine the total number of SN explosions, and the rel- 
ative contributions of different SN types, to the total 
enrichment of the ICM using high-resolution X-ray ob- 



2 



servations. The method relies on using the full set of 
nucleosynthetic yields produced by recent SN models for 
elements with emission lines in the X-ray bandpass to 
probe the total number of SN explosions and relative 
contribution of SNe la. We validate this approach using 
CCD resolution XMM-Newton European Photon Imag- 
ing Camera ( EPIC) and high-resolution RGS observa- 
tions of the cooling flow cluster A3112. A brief summary 
of the method is given in Section [2j The application of 
the method to analysis of XMM-Newton RGS, EPIC, 
and simulated Astro-H Soft X-Ray Spectrometer (SXS) 
data is outlined in Section [3} W e provide our discussions 
and conclusions in Sections |4j and [5] 

2. XSPEC SUPERNOVA ABUNDANCE MODEL (snapec) 
2.1. A Brief Description of snapec 

Most of the metals from oxygen through nickel residing 
in the ICM were synthesized by SN la or SN cc explo- 
sions. We provide a new XSPEC model called snapec 
to determine the total number of SN explosions and rel- 
ative contributions of SN types to the metal content of 
the ICM using elemental yields provided by SN nucle- 
osynthesis calculations in the literature. The expected 
elemental yields of SNe type la, and SNe cc from massive 
(10—50 M ) progenitor stars of various metallicity (0 — 
1 A©) have been intensively studied by severa l authors 
(|Woosley fc Weaver||1995 Iwamoto et al.||1999 Nomoto 
et al.|2 006). For instance, SIN models predict that SJNe la 
produce significant amounts of iron (Fe) , nickel (Ni) and 
silicon (Si) while, for SNe cc, large quantities of oxygen 
(O), neon (Ne) and magnesium (Mg) are produced but 
very little Fe and Ni escapes the compact remnant. 

snapec calculates the mass of the «th element (Mf Ne ) 
in terms of the number of SNe la (N Ia ) and SNe cc (N ) 
explosions that enrich the ICM and the yields per SNe 
la (y! a ) and SNe cc (yf ), 



M. 



SNe 



= N SNe (l 



■R)- x \Ryl a + 



(i) 

(2) 



where R is the ratio of SNe la to SNe cc (R = N Ia /N cc ), 
and N SNe is the total number of SN explosions N SNe = 
N Ia + N cc . The SN yields y\ a and yf are obtained from 
published SN models and stored in the snapec database 
in solar units. The database includes yields of thirty ele- 
ments from SN la explosions of slow deflagration (W7, C- 
DEF), delayed detonation (WDD ; CDDT, ODDT) and 
core 




50 Mo 



an extensive range of progenitor masses (10 
metallicities (0—1 ti mes solar) flTsujimoto et al. 1995 
Iwamoto et "alj[l999l |Woosley fc Weaver||1995| |lNomoto 



et al.||2006[ ) XWW95 and N06, hereafter), snapec adopts 
solar abundances for metals with atomic number < 7 
(e.g., helium, lithium, beryllium, boron, carbon, and ni- 
trogen). 

The snapec model allows a selection of solar abundance 
standards — the datab a se incorporates the sets of | An- 



der s i fc Grevessel ( |1989[ ), |Loddersj ( f2003[ ) , or | Asplund et 
'( 2009p . Since the majority of hydrogen and helium 



al 



in the ICM was produced during the big bang, the stan- 
dard hot big bang nucleosynthesis (BBN) model predicts 
the primordial heliu m abundance by unit mass Yp to 
be 0.2565 ± 0.0010 ( |Izotov fc Thuan||2010| which cor- 
responds to 8.3% helium atoms m number assuming .3 
sola r abundances for met als (e.g., |Markevitch| (2007)). 
The Asplund et al. ( 2009 ) solar abundances we use here 
predict 8.5% of helium which is consistent with BBN pre- 
dictions. 

The SNe cc yields of the ith element, averaged over the 
LMF(</>), is 



< vT >= 



$21 yrOMra) dm 
XT, 2 <M m ) dm 



(3) 



where m\ and mi are the lower and upper limits for the 
masses of SN cc progenitors. Here we assume that stars 
with masses M < 10 M and M > 50 M Q do not explode 
as SNe cc. In this work we use th e standard Salp eter IMF 
slope above 10 M , <f> = m -235 ( |Salpeter||l955[ ). 

It is straightforward to transform Equation (2) to an 
expression for the abundance, by number relative to hy- 
drogen, of the ith element. As long as the total mass 
fraction of metals is small and the He abundance is so- 
lar, 



7,SNe 



z 



M, 



SNe 



N SNe 

Mf(l + i?) 



[Ryl a +<yi c >]• 



(4) 
(5) 



M® is the mass of the zth element for an ICM with solar 
abundances, i.e. 



M. 



© 



M ICM Zf , A z 



m N 



(0 



where Z!f is the abundance of the ith element in so- 
lar units, Ai is the atomic mass of the ith element, and 
mjv = 2.27 x 10 -24 g /H-atom is the averag e nucleon mass 



per hydrogen atom (Asplund et al. 2009). 



This XSPEC model, snapec, has five tree parameters: 
gas temperature (fcT e ), the total number of SN explo- 
sions (N SNe ) and the SN ratio (R) that enriches the 
ICM in the spectral region being analyzed, the cluster's 
redshift, and the XSPEC normalization factor identical 
to those defined for other thermal plasma models (see be- 
low). In addition, there are two fixed parameters, numer- 
ical indices that specify the SN la and SN cc models from 
which the yields are drawn. The predicted spectrum is 
constr ucted using the sin gle temperature thermal model 
(apec) ( |Smith et al~1|2001[ ) and AtomDB 2.0.1 (A. Foster 
et al. 2012, in preparation) to determine the continuum 
and line emission for the metal abundances calculated 
using Equation (5). This method has the advantage of 
simultaneously providing a single set of uncertainties on 
the SN ratio and the total number of SN explosions that 
enrich the ICM since the formation of the cluster. The 
results presented in this work are obtained using the |As-| 
plund et al. (2009) solar abundances. 

snapec redefines the quantities in Equation (5) to ex- 
press the abundances in terms the total number of SNe 
per 10 12 Af Q of ICM plasma — i.e., rescaled to yield val- 
ues appropriate for cluster cores. The total number of 
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Fig. 1. — First— ord er and second— order RGS sources and background spectral fits using the snapec model. The fit was performed using 
|Iwamoto et al.| | |1999[ l SN la yields obtained from W7 models and SN cc yields obtained from the weighted average over a Salpeter IMF. 
The best— fit temperature, total number of supernovae, and SNe type la percentage of the total SNe are 3.69±0.23 keV, (1.24±0.33) X 10 9 , 
and 34.3%±6.4%, respectively (Table [21. The total x 2 is 2818.7 for 2615 degrees of freedom for the best fit to the source spectra, and 
2026.4 for 1498 degrees of freedom for The background fit (Table[2]|. 



SNe is derived by multiplying by the ICM mass (in units 
of 10 12 M Q ) within the spectral extraction region. This 
may be estimated from the XSPEC normalization, which 
is proportional to the emission measure, 



N = 



io- 



n e nndV, 



(7) 



4.7T (1 + Z) 

where n e is the electron number density, Da is the an- 
gular diameter distance, and z is the cluster's redshift. 



TABLE 1 

XMM-Newton RGS Observations of A3 112 



Obs Id 


R.A. 


DEC 


Exposure 


Clean 








Time 


Time 




(J2000) 


(J2000) 


(ksec) 


(ksec) 


0603050101 


03 17 57.4 


- 44 14 14.7 


118.9 


100.1 


0603050201 


03 17 57.4 


- 44 14 12.8 


80.8 


80.4 



3. APPLICATIONS TO X-RAY SPECTRA 

Accurately determining the elemental abundances is 
crucial for understanding the relative contributions of 
different SN types to the metal enrichment of galaxies 
and intra-cluster gas. Earlier attempts to determine the 
relative fraction of SNe la using XMM-Newton measure- 
ments of individual elements such as Si, S, Ar, Ca, Fe, 



and Ni in the ICM \ Werner et al.| 2006b[ |de Plaa et al. 
20071) concluded that T95 and 199 SNe models overes- 



timated Ar and Ca abundances, and cou ld not produce 
accept able fits to XMM-Newton data. |de Plaa et al.| 



( 2007 1 also inferred a number ratio of SNe cc to SNe la 
(using W7 yields) of ~3.5 based on XMM-New ton ob- 
servations of 22 clusters. Grange et al. (20111 studied 
EPIC and RGS observations of two groups of galaxies, 
NGC 5044 and NGC 5813 and found that -30%-40% 
of SNe contributing to metal enrichment of the intra- 
galactic gas were type la while the other ~60%— 70% 
of supernovae then were core collapse. This result was 
consistent with what was previously found fo r M87 by 



Werner et al. (2006a). Simionescu et al. (20091) reported 
that WDU3 or W7 SN la models are needed in order to 
reproduce the low Si abundance in Hydra A clusters of 
galaxies. They also noted that abundance patterns in 
the XMM-JVewton observation of all other clusters pre- 
ferred WDD1 and WDD2 models. However these studies 
proceeded by first determining the best-fit values and un- 
certainties on individual elemental abundances, followed 
separately by fitting abundance ratios to determine the 
relative fractions of SNe la and SNe cc using the T95 
and 199 models. By contrast, in this work we apply the 
yields obtained from SN models to immediately predict 
the spectrum, skipping the intermediate step required 
in previous studies. The best-fit R and N SNe parame- 
ters, and parameter confidence levels, are thus directly 
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TABLE 2 

The Best-fit Parameters of the snapec Model Obtained Using T95, 199, and M10 SN Yields and XMM-Newton RGS Spectra. 



SN la 


SN cc 


kT e 


N SNe 


R 


-ttty 


\ 2 of Source 


X 2, of Bkgd 


Model 


Model 


(keV) 


(xlO 9 ) 




(%) 


Spectra 
(2608 dof) 


Spectra 
(1490 dof) 



T95 (W7) 


T95 


3 


69 


± 





23 


1 


06 


± 





33 





52 


± 





09 


34.3 ± 6.4 


2817 


8 


2026 


8 


199 (W7) 


199 


3 


70 


± 





23 


1 


06 


± 





34 





52 


± 





09 


34.2 ± 6.3 


2818 


7 


2026 


4 


199 (W70) 


199 


3 


74 


± 





23 


1 


07 


± 





33 





50 


± 





09 


33.5 ± 6.2 


2822 


5 


2025 


5 


199 (WDD1) 


199 


3 


81 


± 





23 


1 


15 


± 





35 





57 


± 





10 


36.2 ± 6.8 


2835 


8 


2022 


9 


199 (WDD2) 


199 


3 


80 


± 





22 


1 


09 


± 





34 





48 


± 





08 


32.3 ± 5.9 


2831 


6 


2023 


6 


199 (WDD3) 


199 


3 


80 


± 





21 


1 


06 


± 





33 





44 


± 





08 


30.3 ± 5.4 


2829 


5 


2024 


1 


199 (CDD1) 


199 


3 


81 


± 





22 


1 


17 


± 





35 





59 


± 





12 


37.1 ± 7.1 


2836 


6 


2022 


8 


199 (CDD2) 


199 


3 


81 


± 





22 


1 


08 


± 





34 





45 


± 





08 


31.2 ± 5.7 


2831 


5 


2023 


7 


M10 (W7) 


199 


3 


70 


± 





23 


1 


06 


± 





33 





51 


± 





09 


33.9 ± 6.4 


2819 


2 


2826 


2 


M10 (CDEF) 


199 


3 


66 


± 





30 


1 


28 


± 





43 


1 


23 


± 





27 


56.3 ± 13.7 


2809 


5 


2026 


9 


M10 (CDDT) 


199 


3 


33 


± 





23 


1 


17 


± 





41 


1 


34 


± 





29 


57.2 ± 14.5 


2806 


4 


2026 


8 


M10 (ODDT) 


199 


3 


71 


± 





23 


1 


12 


± 





34 





61 


± 





11 


37.8 ± 7.2 


2823 


6 


2025 


1 



estimated from X-ray spectra and their statistical uncer- 
tainties for a given set of SNe yields, allowing one to use 
the goodness-of-fit to assess the accuracy of the predic- 
tions and compare different yield sets on a well-defined 
statistical basis. 

3.1. XMM-Newton RGS Observations of A3112 

XMM-Newton^ high-spectral resolution RGS instru- 
ment (wavelength bandpass ranging from 5 to 38 A) 
can resolve strong X-ray lines and measure elemental 
abundances accurately in the compact central regions 
of nearby bright clusters of galaxies. We examined the 
RGS observations of the cool core cluster A3112 (see Ta- 
ble [T]) using version 11.0 of the Science Analysis Sub- 
system (SAS) software. For details about RGS data 



reduction and analysis see Bulbul et al. (2012, hereafter, 
Paper I) . We fit the combined RGS first— order and 
second— order spectra using an absorbed snapec model 
in the 7 — 22A and 7 — 16A wavelength range, re- 
spectively. The absorption column density was fixed to 
the Leiden/ Argentine/Bonn Galactic value (Kalberla et 
ail|2005|). The N SNe , R, redshift, kT e , and normaliza- 



tion parameters were allowed to freely vary. The best- 
fit model parameters obtained from the snapec fits to 
RGS spectra for various SN models are shown in Table 
[2l Since, as in Paper I, fitting is conducted using the 
(J-statistic and the C-statistic does not provide a direct 
estimate of the goodness-of-fit, we asses the goodness- 
of-the fit using the corresponding x 2 value. These final 
X 2 values, obtained from the best fit determined by C- 
statistics, are reported in Table [2] 

We calculated the projected gas mass of A3112 con- 
fined in the central 38" region corresponding to the ef- 
fective RGS aperture (Paper I) based on the best-fit nor- 
malization parameter in the single— temperature thermal 
model (apec) fit to the RGS spectra using Equation 
The apec normalization of 1.28 x 10 -2 gives a projected 
gas mass of (4.80 ± 0.70) xlO 11 M within the 38" (~ 52 
kpc) region. We also calculated a spherical gas mass of 
(4.06 ± 0.43) x 10 11 M Q within the 38" (~ 52 kpc) which 
is consistent with the projected mass obtained from the 
apec normalization. This spherical gas mass was calcu- 
lated using the deprojection analysis and modeling of the 
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Fig. 2. — Close-up view of the Fe-L bandpass. The fit was 
obtained by applying the T95 SN la (W7, slow deflagration) and 
SN cc (weighted over 10 — 50 Mq using the Salpeter IMF function) 
yields to the XMM-Newton RGS spectra. The best-fit model 
parameters and the goodness of fit are shown in Table [2] 



XMM-Newton E PIC observations of A3112 (see |Bulbul 
et al. (2010 2012 1). Thus, a conversion factor 0.41 is ap- 
plied to the SNe per 1O 12 M0 derived from spectral fits 
to obtain the total number of SNe, N SNe parameters, 
which are also reported in Table [2] 

We also convert the R parameter obtained from snapec 
fits into the fractional contribution of SNe la to the total 
number of SN explosions as follows: 



Tfla 

H % - 



SNe type la 
(SNe type la + SNe cc) 

R 
1 + R' 



(8) 
(9) 



We find that 30.3%±5.4% - 37.1%±7. % of the total 
SN explosions are SNe la for models that use W7, CDD, 
and WDD yields. Application of M10 CDDT and CDEF 
SN la models implies a significantly higher fraction of SN 
la explosions (~ 57%±14.5%). The reason for this differ- 
ence is that M10 CDDT, CDEF SN la models have lower 
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TABLE 3 

The Best-fit Parameters of the snapec Model Obtained Using T95, 199, and M10 SN Yields and XMM-Newton EPIC Spectra. 



SN la 


SN cc 


kT e 


N SNe 






it 




n % 




X 


ivioaei 


ivioaei 


(Kev ) 


(X10 , 










CV \ 
(/o) 




(^luoo dot J 


1 yo [ W ( ) 


i yo 


q a n _l n no 
0.49 ± U.Uz 


J..UD ± U 


no 


n a n 
U.4U 


i 

± 


n m 


no n i 
ZS.O ± 


i n 
l.U 


Zio/.l 


Tan /"\X?" f 7\ 

iyy { w * j 


TOO 

199 


q Kn _i_ n no 
O.OU ± U.Uz 


l.UO ± U 




r\ Af\ 
U.4U 


i 

± 


n m 
U.U1 


no i 


i n 
l.U 


07Q7 Q 

z < o / .0 


TQQ ?W7D^ 


199 


q m _|- n no 


1.06 ± 


02 


0.39 




0.01 


28.1 ± 


1.0 


2723.9 


199 (WDD1) 


199 


3.47 ± 0.02 


0.82 ± 


02 


0.76 


± 


0.03 


43.2 ± 


2.9 


2893.6 


199 (WDD2) 


199 


3.49 ± 0.02 


0.89 ± 


02 


0.51 


± 


0.01 


33.8 ± 


1.0 


2783.9 


199 (WDD3) 


199 


3.47 ± 0.02 


0.91 ± 


11 


0.42 


± 


0.02 


29.6 ± 


1.9 


2749.8 


199 (CDD1) 


199 


3.41 ± 0.01 


0.82 ± 


11 


0.81 


± 


0.03 


44.8 ± 


2.9 


2927.6 


199 (CDD2) 


199 


3.50 ± 0.02 


0.87 ± 


11 


0.49 


± 


0.02 


32.9 ± 


1.9 


2767.1 


M10 (W7) 


199 


3.51 ± 0.02 


1.06 ± 


01 


0.40 


± 


0.01 


28.6 ± 


1.0 


2744.1 


M10 (CDEF) 


199 


3.54 ± 0.02 


1.36 ± 


11 


0.77 


± 


0.03 


43.5 ± 


2.9 


2775.8 


M10 (CDDT) 


199 


3.39 ± 0.02 


1.11 ± 


10 


1.28 


± 


0.07 


56.1 ± 


6.5 


3627.5 


M10 (ODDT) 


199 


3.49 ± 0.02 


0.83 ± 


13 


0.82 


± 


0.04 


45.1 ± 


3.8 


2842.3 



Fe abundance compared to the W7, W70, CDD, WDD, 
and ODDT models. Since Fe is mainly produced by SN 
la explosions, a lower abundance of Fe in SN la models 
in M10 CDDT, CDEF yields results in a higher number 
of SN la explosions and thus higher SN la fraction. The 
fraction of SN la explosions inferred using the 199 CDD 
and WD D models reported in Tableware consistent with 
those in de Plaa et al. (2007) based on elements acces- 
sible to the FPlC~detectors : Si, S, Ar, Ca, Fe, and Ni. 
However de Plaa et al. ( 2007[ ) find a lower SNe la fraction 
for W7 SN la models. 

The SN models we use in this work produce equally 
good fits to RGS data (see Table [2]). Therefore, global 
fits to the currently available high— resolution X-ray ob- 
servations of clusters of galaxies do not allow one to 
distinguish between different types of SN models in an 
individual cluster. The RGS first— and second— order 
spectral fits obtained using the 199 SN la (W7) yields 
and Salpeter-IMF-averaged 199 SN cc yields are shown 
in Figure [T] We also show a close-up view of the fit in the 
RGS Fe-L bandpass obtained using the T95 SNe la (W7) 
and Salpeter-IMF-averaged T95 SN cc yields in Figure 
[2j Figure [2] shows that the T95 SN yields accurately 
reproduce the RGS spectrum in the Fe-L bandpass. 

3.2. XMM-Newton EPIC Observations of A31 12 

We examined the EPIC observations of the cool core 
cluster A3112 (see Table 1) using version 11.0 of the SAS 
software. The EPIC data processing and background 
modeling were carried out with the XMM-Newton Ex- 
tended S ource Analysis Software (XMM-ESAS) and 
methods flKuntz fe Sn owden|[2008) lif" 
since the EF 

strain Si, S, Fe and O abundances only in the innermost 
30" region (see Paper I), we examined only the central 
30" spectra in this work. We fit the cluster emission using 
the absorbed snapec model allowing the absorption col- 
umn density, N SNe , R, redshift, fcT e , and normalization 
parameters freely vary. The best fit parameters of the 
instrumental and cosmic X-ray background model com- 
ponents were fixed to the best-fit values obtained from 
the overall fits and reported in Paper I. A conversion 
factor 0.41 was again applied to the SNe per 10 12 A/ Q 
derived from EPIC spectral fits to obtain the N SNe pa,- 



( |Kuntz fe Snowdeli||2008| ISnc-wden et aL][2008 | 
described in Paper 1. Since the EPIC spectra can con- 




Energy (keV) 



Fig. 3.— snapec fit to the XMM-Newton MOS1 (red), MOS2 
(blue), and PN (magenta) spectra from the innermost 30" region 
of A3112 obtained from the observation 0603050101. The back- 
qr ound model paramete rs were fixed to the best fit val ues reported 



by\Bulbul et al. \2012\. The fit was performed usi ng \Iwamoto et\ 
[at| \1999\ WW SN la and \Iwamoto et al.\ \1999\ SN cc yields 
obtained from the weighted average yields over a Salpeter IMF nu- 
cleosynthesis yields. The best fit temperature and SN ratio (R) 
and the total number of SN (N SNe ) parameters, 3.50 ± 0.02 keV, 
0.39±0.01, and (1.06±0.02) X 10 9 , respectively, are shown in Table 




rameter. The best-fit model parameters obtained from 
the snapec fits to EPIC spectra for various SN nucle- 
osynthesis models are shown in Table [3| These final % 2 
values, reported in Table [3j were obtained from the best 
fit determined by C-statistics. 

The temperature, SN ratio (R), and the total num- 
ber of SN measurements obtained from EPIC spectra 
are consistent with RGS results (Table [2]) at the ler 
level. The la difference is due to different temperature 
measurements reported by EPIC and RGS observa- 
tions. Most of the SN nucleosynthesis models we use in 
this work produce equally good fits to EPIC data as 
in the case of RGS analysis. High signal-to-noise CCD 
resolution EPIC spectra yield stringent constraints on 
the model parameters (e.g. SN ratio (R) and the total 
number of SNe) but cannot distingui sh between different 
types of SN nucleosynthesis models. Maeda et al. (20101 



G 




Fig. 4. — Simulated Astro-H SXS spectra using the snapec model. The fit was performed using |Twamoto et al.| ( |~1999[ ) SN la yields 
obtained from W7 models and SN cc yields obtained from the average yield computed using Salpeter IJVLb 1 weighting (as in[TJ|. The red 
color shows the best-fit snapec model with T95 W7 SN la and SN cc yields. The contribution of SNe la explosions to the metal abundances 
are shown in green while the contribution of SN cc explosions are shown in blue. 



CDDT SN la m odels yield a worse fit t o the EPIC spec 



tra compared to Iwamoto et al. (19991; Tsujimoto et al 



p95 ) or |Maeda et al.|fl20l0| ) W77TIDEF, ODDT SlN la 
models. 

We find that fits to the EPIC spectra find a slightly 
wider range of SN la fraction, so that 28.1%±1.0% — 
44.8%±2.9% of the total SN explosions are SNe la for 
models that use W7, CDD, and WDD nucleosynthesis 
yields. As in the case of the RGS analysis, the M10 
CDEF, CDDT SNe la models imply a significantly higher 
fraction of SN la explosions (43.5±12.9%-56.1±6.5%) 
as a result of lower Fe abundance compared to the W7, 
W70, CDD, WDD models. We also found that all W7 
models from several simulations (e.g. T95, 199 and M10) 
produce similar goodness-of-the-fit, the total number of 
SNe, and SN la fraction. The snapec fits to EPI C MOS1 
MOS2 , and PN spectra performed using the |Iwamoto] 
et al.| dl999]) W70 SN la and Salpeter-IMF-weighted 



Iwamoto et al. (1999) SN cc metal yields are shown in 
Figure |5J 

3.3. Astro-H Simulations of A3112 

The Japan/U.S. Astro-H Observatory is scheduled to 
launch in 2014 and will carry an X-ray calorimet er, the 
SXS, with a high energy resolution of ^4.5 eV (|Taka 



|hashi et al.||201f}| ). The detector will be able to deter 
mine line widths and metal abundances with high preci- 
sion, revolutionizing our understanding of physical pro- 
cesses (e.g. turbulence, the contributions of SN explo- 



sions to metal enrichment) in the ICM. We simulated 200 
ks SXS observations of A3112 using the X-ray events sim- 
ulator software (simx 1.2.1), based on the temperature 
and abundance measurements obtained from the snapec 
fits convolved with 199 SN la (W7) and 199 Salpeter- 
IMF-weighted SN cc yields, simx uses predefined detec- 
tor response matrix, effective area, and predicted back- 
ground for the SXS instrument to produce an event file 
by convolution with the telescope's point spread func- 
tion. 1 The SXS spectrum is extracted using the FTOOL 
xselect based on the event file produced with the simx 
software. 

We then fit the simulated SXS spectrum with the 
snapec model convolved with various SN la and SN cc 
yields for variable temperature, N SNe , R, redshift, and 
normalization. Figures [4] and [5] show the soft X-ray band 
and a close-up view of the Fe-K bandpass for fits to the 
simulated SXS spectrum of A3112 obtained assuming the 
same SN yields used to create the spectrum. The relative 
contributions of SN la and SN cc products to the spec- 
trum are shown as green and blue lines in Figures [4] and 
[5j respectively. The figures illustrate how the significant 
amounts of Fe and Mg produced by SN la explosions, and 
large quantities of O synthesized in SN cc explosions, are 
manifest in the spectrum. The best-fit parameters and 
goodness-of-fits obtained using the 199, T95, and M10 
SNe yields are shown in Table [4] The features from met- 

1 http:/ /hea-www. harvard.edu/simx/index. html 
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200 ks Astro-H 
simulations of A3112 
Fe-K bandpass 



SNe la 

SNe cc 
SNe la+SNe cc 




6.2 

Energy (keV) 



Fig. 5. — Close-up view of the Fe-K bandpass of the 100 ks simu- 
lated Astro-H SXS spectrum, based on the snapec model convolved 
with 199 SN la and SN cc yields. The red color shows the best-fit 
snapec model for T95 W7 SN la and SN cc yields. The contri- 
bution of SN la explosions to the metal abundances are shown in 
green while the contribution of SN cc explosions are shown in blue. 



als primarily produced by SNe la explosions such as S, 
Ar, Ca, and Ni that were inaccessible to the EPIC and 
RGS now become important diagnostics for Astro-H data 
analysis. As a result of the accurate abundance measure- 
ments of these metals, uncertainties in the N SNe and R 
parameters are reduced, and they may be determined 
more robustly and with higher precision (see Table [4 ) . 

As expected, good fits arc found using SN models that 
have similar elemental yields, e.g. using T95 (W7) or 
199 (W70) SN la models, to those of the 199 (W7) SN 
la and SN cc yields adopted to produce the simulated 
Astro-H spectrum (see Table EJ). On the other hand, the 
significantly higher \ 2 values obtained from fits with SN 
models that have markedly different yield patterns, e.g. 
199 CDD and 199 WDD models, demonstrates that future 
Astro-H observations of galaxy clusters will allow one 
to distinguish between different SN models and provide 
fundamental constraints on SN nucleosynthesis. 

4. DISCUSSION 

By developing and applying a new spectral analysis 
methodology to RGS spectra, we measure the total num- 
ber of SNe that have enriched the ICM presently in the 
inner 52 kpc of A3112 and the relative contributions of 
SNe la and SNe cc. In this section we compare this 
with the enrichment expected based on the number of 
SN la and SN cc explosions derived from the amount of 
starlight and well-established empirically— measured SN 
rates per star formed. We also investigate prospective 
local enrichment sources related to the brightest cluster 
galaxy (BCG). In this way we can interpret our data 
analysis in the context of the origin of the gas and its 
metals. 

4.1. Global Cluster Context for the Number of SNe 

Within its virial radius, a sufficiently massive cluster 
is approximately a closed box. For a global efficiency of 
converting gas to stars £ s / , the total mass in stars formed 
is 



where M\, ary is the total baryon mass. At the present 
time the total mass in stars, whether contained in indi- 
vidual cluster galaxies (including the BCG) or associated 
with intra-cluster light (ICL), is 



M* — Mb ar y e s f (I — i?* 
and the mass in gas is 



(11) 



M gas = M hary - M* (12) 
= M hary [l-e sf {l-R*)}, (13) 

where i?* is the mass return fraction. We can safely 
neglect the distinction between the mass in gas and the 
mass in the ICM for a cluster as rich as A31I2 where the 
gas mass in galaxies is relatively small, and henceforth 
equate M gas with M\cm- The star formation efficiency, 
in terms of the observable M\cm /M*, 



e sf = (I - i?*)"^! + Mjcm/M*)- 



(14) 



The total numbers of SN cc and Type la explosions 
can be expressed as 



M, 



*form 



(15) 



and 



N Ia = r! Ia M* form , (16) 

where r\ cc and r/ Ia are, respectively, the specific numbers 
of SN cc and SN la explosions per star formed. Connect- 
ing rf c and r\ la to the snapec model parameters defined 
in m 



R = 



and 



N SNe = v {l-R*)~ L M* 



(17) 



(18) 



where rj — r\ + r\ . We note that these refer to all of 
the SNe while in Section[2j N SNe and R parameters refer 
only to SNe that enrich the ICM. 

Adopting a "diet Salpete r" IMF that produces rela- 
tively fewer low ma ss stars ( Bell fc de JongpOOl I , Maoz 
& Mannuccrj ( |201l| ) estimated r/ Ja to be ~ 0.002. Adopt- 
ing the mass return fraction and spe cific SNe cc rate for 
the diet Salpeter IMF : R., ~ 0.35 ( |Fardal et ah] [2007 
O'Rourke et al.||201l|), n Ia ~ 0.002, and n 11 ~ 0.008 
( |Maoz fc Mannucci|26li~| Botticella et al.|2012[ ), one pre- 
dicts R ~ 0.25, which corresponds to i?£, a = 0.2, and 



n ■ 



which corresponds to R% 
0.01. For clusters in the A3112 mass range (e.g.. 
3.0 x 10 14 M Q , Pap er I), recent analysis finds. 



2011 Lin et 



ty pically , M*/Mi C m ~ 0.1 dBalogh et al 
al. 2012). Therefore, from Equation (T5), the estimated 
total number of SN explosions per 10 M of ICM is 
~ 1.54 x 10 9 . Multiplying by the correction factor of 
0.41 appropriate to the core of A3112, N SNe becomes 
~ 6.3 x 10 8 . Comparing with our results in Tabled we 
may conclude that the inner 52 kpc of A3112 has Deen 
enriched by more SNe, with a higher percentage of SNe 
la, than is true globally for a typical cluster of its mass. 
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TABLE 4 

snapec Parameters for the Best Fit to the Simulated Astro-H Spectra of A3112 Assuming 199 SN la and SN cc Yields 



Id 


DIN CC 


i-T 

Me 


1 V 


r> 


X 


Model 


Model 




CxlO 9 ! 




f886 dnfl 


T95 (W7) 


T95 


3.77 ± 0.03 


1.10 ± 0.07 


0.50 ± 0.02 


1124.5 


199 (W7) 


199 


3.77 ± 0.03 


1.09 ± 0.06 


0.51 ± 0.02 


1119.9 


199 (W70) 


199 


3.76 ± 0.03 


1.08 ± 0.06 


0.49 ± 0.02 


1113.8 


199 (WDD1) 


199 


3.78 ± 0.03 


1.06 ± 0.06 


0.60 ± 0.03 


1319.4 


199 (CDD1) 


199 


3.69 ± 0.05 


1.07 ± 0.08 


0.33 ± 0.02 


1348.9 


M10 (CDEF) 


199 


3.76 ± 0.03 


1.46 ± 0.08 


0.89 ± 0.05 


1223.9 


M10 (CDDT) 


199 


3.83 ± 0.03 


1.16 ± 0.05 


1.54 ± 0.12 


1795.0 


M10 (ODDT) 


199 


3.84 ± 0.03 


1.07 ± 0.07 


0.69 ± 0.03 


1292.6 



4.2. SN Metals Locked Up in Stars 

The galactic mass in clusters is dominated by early- 
type galaxies that form their stars quickly. This results 
in the well-established enhancement in [a/Fe], the abun- 
dance ratio of a-elements to Fe (expressed as the log- 
arithm with respect to solar) — i.e., SNe cc are pref- 
erentially locked up in stars. In their investigation of 
the giant e lliptical galaxy NGC 4472, |Loewenstein fe| 
Davisl ( |2010[) found that a ratio of SNe la to total SNe 
■^fJjiaVjj^Ne* ^ on ( N ccy N SNe* _ 0.89) and a 

number of supernova per mass in (present-day) stars 
of N SN /M* ~ 0.0083, resulted in [a/Fe]* - 0.25 and 
Zp ele ~ 1 (as observed in this particula r galaxy, but typ- 



ical of the class) assuming yields from Kobayashi et al. 
(2006). This enables us to estimate the lock-up correc- 
tions, r] cc * and n Ia * , to the total specific numbers of SNe 
per star formed. This provides, in turn, the correspond- 
ing values available to enrich the ICM, t)\q M and nf CM : 



V = 



= 6.0 x \{r 4 Z Fe ^ 



N c 



A'L 



:4.8 x 10~ 3 Z Fe>k , 



(19) 
(20) 

(21) 
(22) 



Vwm = V Ia - V Ia * = 1-4 x 10- 3 Mq\ (23) 



and 



vTcm=V cc -V cc * = 3.2x10-^. 



(24) 



That is, the metal production from ~ 60% of SNe cc, 
and ~ 30% of SNe la, must be locked up in stars to en- 
rich them to solar Fe abundances and [a/Fe]* ~ 0.25. 
This results in revised predictions of R ~ 0.30 — consis- 
tent with our RGS results for A3112, but an even lower 
expected total number of SNe - N SNe ~ 2.9 x 10 8 SNe. 

4.3. Direct SNe Injection 



Sand et al. (2012) recently estimated a specific SN la 
rate in cluster galaxies of rsNia ~ 0.04 SNuM 2 within 



-R200 (~ 1 Mpc) (see also Maoz et al. (2010)) that may 



be directly injected into the ICM by the r < 52 kpc 

1 SNuM = Supernova rate per 100 yr per 1O 1O M0 in stars 



stellar population. The implied level of SNIa available 
for enrichment accumulated over time r is 



N 



= 10~ 4 AF(< 52 



kpc) ^ 



rSNIa 



)(: 



(1.1 S.V//.W I \ WPyr) 



(25) 

where the stellar mass refers to the r < 52 kpc region, 
and rsNia should now be interpreted as an average over 
r. We use the Hernquist approximation (Hernquist 1990) 
to deVaucouleurs profiles for the A3112 BCG and 1CL 
components and estimate = 3.6 x 10 11 M Q within 52 
kpc — i.e., M* = 0.88MicAf. A large star-to-gas ratio, 
M* ~ M\cm an d a long timescale r > 10 9 yr (with the 
accompanying increase in r$Nia) would be required for 
this to be a significant source of enrichment. 

4.4. The Cluster Enrichment Puzzle as Realized in the 
Core of A3 112 

For most of the yield sets we consider, the level of 
ICM enrichment for the central (4.06 ± 0.70) xlO 11 M 
requires ~ 1.1 x 10 9 SNe, with SNe la accounting for 
~ one-third of the total. Although our SN parameters 
are derived for the inner 52 kpc, the high gas mass and 
temperature imply that this gas is primarily intra-clustcr 
in origin. If we account for the fraction of SN products 
locked up in cluster galaxy stars, which is higher for SNe 
cc based on the measured stellar [a/Fe] enhancement, we 
can explain the proportions of SNe la and SNe cc that 
we infer enriched the ICM. However, the total number of 
SNe we infer is greater than expected by a factor of ~ 4 
for a typical star-to-gas ratio, M^jM\cu ~ 0.1. These 
conclusions are altered by direct injection of SN la ejecta 
in the BCG only if accumulated over a large fraction of 
a Hubble time; and, are exacerbated for the M10 CDEF 
and CDDT models. Evidently, the stellar population 
responsible for enriching this material is characterized by 
enhanced efficiency in producing SNe, and/or the mass 
in stars responsible for the enrichment, M* > 0.3Mtcm- 

Larger estimates of M*/Mic m ~ 0-2 at the A3112 
mass range ma y be fou nd in Gonzalez, Zaritsky, & 



Zabludoff (|2007[) (see, also|Langana et al.|Q201ip). Indeed" 



A3 112 wa s investigated by |Gonzalez, Zaritsky, fe Zablud-_ 
off| ( |2007l ), from which we may derive M*/Mi CM ~ 0.27.' 
The implication is that there is substantial I CL compo- 
nent unaccounted for in other studies (see also |Bregman,| 



Anderson, & Dai (20101)) and that the star formation 



efficiency is very high Equation (14) 
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Even if of intra-cluster origin, the enrichment of this 
inner cluster gas is likely to reflect the somewhat special 
conditions in the cluster core where high initial ovcrden- 
sity may result in enhanced efficiency of galaxy and star 
formation. Thus, the effective value of M^/Micm m &y 
very well be larger than is representative of the ICM as 
a whole. At the same time, for rich galaxy clusters in 
general, the global ICM shows enrichment beyond what 
is ex pected based on the stars we see in ga l axies to- 



al. 



day (Portinari et al.||2004[ |Loewenstein||2006 Maoz et 



2010; BrcgmanT~a.nderson, &: Dai||2010|~ 



Another possible mechanism contributing to the 
greater than expected number of SNe found in our fits 
is strong early enrichment of the lowest-entropy X-ray 
emitting gas stripped from early galaxies. The heavier 
metals found in gas stripped at large radii will settle into 
the cluster core due to the central gravitational poten- 
tial, while the stars that produced them remain at larger 
radii. Therefore, a significant fraction of metals residing 
in the cluster core that we observe today ma y not have 
ori ginate d from the local stellar population (iMillion et 
al. 2011[ ). In this scenario, the apparent total expected 
number of SNe will be enhanced relative to that inferred 
from / M icm ■ A more extensive study of a larger 
sample is needed to fully address this issue. 

5. CONCLUSIONS 

In this work we introduce a new approach to deter- 
mine the total number of SNe explosions and the rela- 
tive contribution of SNe la to the total enrichment of the 
ICM using high— resolution X-ray observations of clus- 
ters of galaxies. The method has been implemented as 
an XSPEC model called snapec. The method relies upon 
the apec model for spectral emission, but uses the metal 
abundances determined from SNe yields from the latest 
published SN la and SN cc e xplosion models normalized 
to the so lar abundances (e.g.JAnders &: Grevesse (1989), 



Lodders| ( |2003 
turns the tota 



Asplund et al.|(|2009p ). This modelre 



number of SJN explosions and the relative 



contribution of the SN types in the spectral extraction 
region that enrich the ICM integrated over the cluster 
life-time, although it does not explicitly constrain the SN 
products locked up in stars in member galaxies. Contrary 
to previous studies in the literature, this method has the 
ability to use all relevant elements, allowing spectra and 
SN yields obtained from published models to determine 
the best fit to their abundances. The main advantage of 
this approach is that it provides a self-consistent single 
set of uncertainties on the total number of SN explo- 
sions and relative fraction of SN types by directly al- 
lowing these parameters to be determined by SN yields 
provided by simulations. This method provides both rig- 
orous assessments of uncertainty and improved statistical 
efficiency in estimating the total number of SN explosions 
and relative SN fractions. 

We use XMM-Newton EPIC and RGS observations 
of A3112 to validate the method and determine the total 



number of SNe a nd fraction of SNe la us i ng various pub- 
lished SN yields (Tsujimoto et al. 1995 Iwamoto et al. 



1999 Mae da et al.||2010| ). Since the T95, 199, and M10 
SJN models produce equally good fits to RGS data, cur- 
rent CCD resolution EPIC and high— resolution RGS 
observations of clusters of galaxies do not allow one to 
distinguish between different SN models. The SN la frac- 
tion reported by EPIC observations are consistent with 
RGS measurements at la level. We also find that T95, 
199, and M10 models produce consistent values of the 
total number of SNe and SN la fractions. 30.3% ± 5.4% 
to 37.1% ± 7.1% of the total SN explosions are SNe la 
and the total number of SN explosions are in the range 
of (1.06 ± 0.34) x 10 9 to (1.28 ± 0.43) x 10 9 . The best-fit 
SN la ratios for WDD and CDD we obtained from RGS 



spectra are consisten t with |de Plaa et al 



( |2007| ) mea- 
find a lower 



surements. However, de Plaa et al. 
SN la fraction for W7 SJN la models'. 

We also applied this method to simulated Astro-H SXS 
observations of A3112. Astro-H simulations were per- 
formed using 199 W7 SN la and 199 Salpeter— weighted 
SN cc models. Assuming that 199 W7 models reflects 
the true nature of the ICM enrichment, higher \ 2 values 
obtained from the 199 WDD, 199 CDD, and M10 CDDT 
models indicate that we will be able to distinguish be- 
tween different SN models with Astro-H observations of 
clusters of galaxies. 

The estimated mix of SNe la and SNe cc (30% and 
70%, respectively) expected in the ICM of A3112 based 
on well-established, empirically measured SN rates per 
star formed is close to what we infer to have enriched 
the ICM in the inner 52 kpc of A3112 based on this new 
method of high resolution spectral analysis. The inferred 
level of SN enrichment, - 2.7 x 10 9 SNe per lO^M" 1 of 
ICM, corresponds to M*/Mjcm ~ 0.3-0.4 — compared 
to M*/M\cm ~ 0.1 that has been estimated as typical 
for clusters in the A3112 mass range. 

When applied to RGS observations of a large sample of 
cool core galaxy clusters, this method can provide strin- 
gent constraints on the contribution of SN explosions to 
the metal enrichment of the ICM and the evolution of the 
total number of SNe la with redshift. This method will 
also provide an independent measure of a key diagnostic 
of the nature of SN la progenitors, the dependence of the 
evolution of SN la rate to the distribution of the delay 
time between the formation of the progenitor system and 
its explosion as an SN. 
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